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Abstract Phylogeographical analysis was conducted on butterflies of the genus Erebia in Bulgaria, which occupies a 
major position on the Balkan Peninsula. We studied the haplotypes of base sequences of 942 bp in total using ND5 (432 
bp) and CO7 (510 bp) which are often used for analysis of intraspecific molecular variation. Samples are 65 individuals 
of 10 Erebia species from northern Bulgaria (Western and Central Balkan Mts) and southern Bulgaria (Rodopi, Pirin and 
Rila Mts), and also 11 individuals of 10 species from central Europe mainly Switzerland and central Asia. 

A geographical distance-based isolation mechanism may be functioning for E. aethiops, E. euryale, and E. medusa, which 
are widely distributed in Europe, with just a one-base substitution in the sequence between the north and south populations 
of E. aethiops, and the other north-south species being in an equivalent situation. On the other hand extremely large 
genetic differences between the north and south Bulgarian populations of E. ligea and E. oeme have been discovered. 
This is probably a result of reproductive isolation between the north and south populations prior to the last glacial period 
which remains till now. Especially with E. ligea, it may be that two lineages, the European and Asian, are facing each 
other in the central Balkan Peninsula. In addition, in the E. alberganus and E. pandrose populations, reproductive isolation 
between the Balkan populations and central European populations seems to have continued since the age prior to the last 
glacial stage, since the genetic differences between them are very large. This suggests that some northern butterflies like 
these two Erebia species, now distributed very sporadically, might not have expanded their distribution and neither mixed 
each other in the last glacial period even on the continent. 
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Introduction on northern organisms, reported cases show more varied 
patterns in terms of such factors as the locations of refugia 
during the warm periods and dispersal routes during the 
ice age than is the case with southern organisms (Hewitt, 


1996; Schmitt, 2009). 


The Balkan Peninsula have been known as a hotspot for 
European biodiversity (Griffiths et al., 2004) since, together 
with the Iberian and Italian Peninsulas, it functioned as a 
refugium for southern (Mediterranean) organisms during 
the ice age (Hewitt, 2000). The Balkan Peninsula generally 
indicates an area south of the Danube river. On the western 
side of the Peninsula, the Dinaric Alps stretch from 
northwest to southeast, and are among the five largest 
mountain ranges in Europe, with Prokletije (2,692 m) as 
the highest peak. At the southernmost end, the mountains 
are divided to two mountain ranges, the Balkan Mountains 
(highest peak Betov, 2,376 m) running from west to east, Balkan Mts 
and the Pindos Mountains (highest peak Katafidi, 2,393 Q 
m) running from north to south. There is a high level of m n 

biodiversity compared to their rather narrow regional areas, 
such as the butterfly fauna with 214 species in Bulgaria 
(Abadjiev, 2001), and 232 in Greece (Pamperis, 1997). 





Pirin Mts 


Rodopi Mts 
③ 


In studies on northern (non-Mediterranean) organisms, 
however, some cases have been reported of organisms 
having small refugia in the far northern region other than 
the Mediterranean, corresponding to the strength of cold 








tolerance of the organisms (e. g. Steward and Lister, 2001; 
Schmitt et al. 2006). Although there have been few studies 


Fig. 1. Balkan Peninsula. Map showing the approximate locations 
of the sampling sites. 
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(e) (f) 
Fig. 2. Habitats of the genus Erebia 
(a) Rila Mts, distant view, (b) Rila Mts, Mt. Borovets (h=2,550 m). E. pandrose on slope areas with short 
grasses and rocky outcrops. (c) Rila Mts, Grqnchar Chalet-Dzhanka Pass (h=2,500 m). E. orientalis and E. 
oeme on grassy slopes near the timberline. (d) Pirin Mts, Bezbog Chalet (h=2,650 m). Many E. rhodopensis 
on the boundary areas of rocky outcrops and short grass areas. (e) Rodopi Mts, Devin, Trigrad (h=1,600 


m). E. euryale and E. ligea. (f) Western Balkan, Chiproviska Planina (h=1,350 m). E. orientalis, E. alberganus, 
E. ottomana on tall grassy and flowery slopes with sporadic low coniferous trees. 
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Fig. 3. Erebia butterflies in the field of Bulgaria 
(a) Erebia ottomana (Pirin Mts, Gotse Delchev-Sandanski, 21 July 2009) 
(b) Erebia orientalis (Western Balkan, Chiproviska Planina, 14 July 2009) 
(c) Erebia euryale (Pirin Mts, Gotse Delchev-Sandanski, 21 July 2009) 
(d) Erebia ligea (Rodopi Mts, Devin, Trigrad, 18 July 2009) 
(e) Erebia alberganus (Western Balkan, Godech, Petrohan Pass, 13 July 2009) 
(f) Erebia rhodopensis (Rila Mts,Borovets, Moussala, 22 July 2009) 
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Schmitt et al. (2006) studied the genetic structure of 
populations of the alpine butterfly E. epiphron over its 
range of distribution. The strongly hierarchical genetic 
structure among populations geographically distant from 
each other led them to surmise that these populations were 
already reproductively isolated in the last ice age rather 
than that the present distribution of local populations had 
been formed in the short postglacial period. Due to the 
extreme dry climate of the ice age, some alpine butterflies 
which preferred humid environments were incapable of 
wide distribution even in the cooler age, and may have 
been sporadically distributed in regions in which humidity 
was highly preserved microclimatically. The same patterns 
are observed in alpine plants (Schonswetter et al., 2003a, 
2003b). 


Haubrich and Schmitt (2007) examined the genetic structure 
of the alpine-endemic butterflies E. sudetica and E. melampus. 
They detected strong differentiation into three lineages 
with the genetic distances between the two E. melampus 
groups being larger than between each of the two E. 
melampus groups and E. sudetica. These results give 
evidence for the existence of a species complex within 
the E. melampus / sudetica group and indicate a discon- 
tinuous distribution within this group during at least the 
last ice age. One of them, E. sudetica, is found in the 
northern Alps and most probably had its refugium north 
of the glaciated Alps. The western E. melampus group 
might have had a refugium at the margin of the southwestern 
Alps, and the eastern group in the lower altitudes of the 
southwestern and/or eastern Alps. 


Erebia species inhabiting the coniferous forest zone are 
considered not to have inhabited environments lacking 
coniferous forests even in the ice age. Schmitt and Haubrich 
(2008) examined the genetic structure of populations of 
forest-dwelling E. euryale in four different European 
mountain systems (Pyrenees, Alps, Carpathians, and Rila). 
They found four different groups (Pyrenees, western Alps, 
eastern Alps, and southeastern Europe) well supported by 
many analysis methods. The genetic diversity of the 
populations was highest in the southeastern European 
group and decreased stepwise westwards. They assume 
that the most important refugia for the coniferous forest 
biome in Europe was located in southeastern Europe, and 
that a fragment of this habitat type survived along the 
foothills of the Alps and Pyrenees. 


Allozyme analysis of E. medusa, more widely distributed 
in the lower altitude areas than E. euryale over large regions 
of Europe, revealed a significant population differentiation 
(Schmitt and Seitz, 2001). Phylogenetic analysis showed 
a division into four major lineages, the eastern group 
including the Czech Republic, Slovakia and Hungary, the 


central group France and Germany, the lineages being 
delimited from western Hungary, and from northern Italy. 
The genetic difference between the populations is small, 
probably arising from disjunct distribution during the last 
ice age. Genetic lineages evolving during glacial isolation 
frequently come into contact as the result of post glacial 
range expansions. Hybridization often occurs along these 
contact zones. Hybrid populations between the eastern 
and the western genetic lineage of E. medusa were detected 
in the Hercynian mountain area along the Czech-German 
border (Schmitt and Muller, 2007). Northern butterflies 
that had less strong cold tolerance were forced into 
vicariance during the ice age before expanding their 
distribution during the postglacial period. 


Evidence of four different Quaternary glacial stages has 
been found in the north-west Iberian Peninsula. The 
phylogeographic history of Iberian E. triaria was analysed 
(Vila et al., 2005). One of the NW Iberian populations, a 
subspecies endemic to the Xistral mountain range, showed 
a high level of genetic divergence from other populations. 
Isolation after an ancient glacial stage, followed by 
allopatric differentiation, may account for such high 
differentiation. The genetic pattern shown by the rest of 
population samples consisted of two lineages, likely 
reflecting that the ancestors of these populations sought 
refuge in at least two different areas during a subsequent 
glacial age. If an analysis is conducted with many more 
sample points, multiple events may be inferred, including 
range expansion and/or fragmentation, as in the case of 
E. niphonica in Japan (Nakatani et al., 2007). 


In this paper we discuss the phylogeographic history of 
the genus Erebia in Bulgaria, the core of the Balkan 
Peninsula, by comparing the alpine butterflies on the 
Japanese islands. 


Materials and Methods 


The authors collected 65 individuals of 10 species from 
the northern mountain systems (Western and Central 
Balkan Mts) and southern mountain systems (Rodopi, 
Rila and Pirin Mts) in Bulgaria, and 11 individuals of 10 
species from central Europe, mainly Switzerland, and 
central Asia. The bodies of the sample specimens without 
wings were preserved in 99% ethanol and kept at 4°C until 
DNA was extracted. 


For genetic markers, we used a part of the NADH 
dehydrogenase subunit 5 (VD5) and cytochrome oxidase 
subunit I (CO) genes in mitochondrial DNA. DNA was 
extracted from the thoracic flight muscles or legs by 
DNeasy Tissue Extraction Kit (Qiagen, Hilden, Germany). 
PCR was conducted in a 104d reaction system by adding 
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Table 1. Haplotype list. 
Detected locality in Bulgaria; WB: West Balkan. CB: Central Balkan, RO: Rodopi Mts, PI: Pirin Mts, RI: Rila Mts 





























: Accession No. Detected locality in Bulgaria 
Haplotype Sampled locality NDS Col "^ We CB R S Ec RI 
E. aethiops 
AE000 Swiss: Goppenstein, Wallis AB674329  AB674366 1 
AEO10 Bulgaria: Chiproviska Planina, West Balkan AB674330 AB674367 2 O 
AE011 Bulgaria: Rozhen Pass near Smolyan, Rodopi AB674331 AB674368 3 O 
E. alberganus 
AL000 Swiss: Zermatt, Wallis AB674332 AB674369 1 
ALO10 Bulgaria: Chiproviska Planina, West Balkan AB674333 AB674370 2 ONES 
E. euryale 
EUOO0 Swiss: Gross Scheidegg, Bern AB674334 AB674371 ] 
EUOIO Bulgaria: Orelyak, Pirin AB674335 AB67432 5 © O: wx O 
EUO11 Bulgaria: Chiproviska Planina, West Balkan AB674336 AB674373 1 〇 
EUO12 Bulgaria: Betov, Central Balkan AB674337 AB674374 1 Q 
EU013 Bulgaria: Pomporovo Resort, Rodopi AB674338 AB674375 1 〇 
EU014 Bulgaria: Trigrad Gorge, Rodopi AB674339 AB674376 2 O 
E. ligea 
LI000 Swiss: Goppenstein, Wallis AB674340 AB674377 1 
LIOO1 Swiss: Maienfeld, Graubunden AB674341 AB674378 1 
L1100 Mongolia: Bayanchandman, Tov AB674342  AB674379 1 
LIO10 Bulgaria: Chiproviska Planina, West Balkan AB674343  AB674380 2 Q 
LIO11 Bulgaria: Chiproviska Planina, West Balkan AB674344 AB674381 1 O 
LIO12 Bulgaria: Chiproviska Planina, West Balkan AB674345 AB674382 3 Q 
LIO13 Bulgaria: Trigrad Osmanov, Rodopi AB674346 AB674383 5 O 
E. medusa 
ME000 Germany: Baden-Wurttemberg AB674347 AB674384 1 
ME100 Mongolia: Terelji, Tov AB674348 AB674385 1] 
ME010 Bulgaria: Central Balkan, Pirin(Mt. Bezbog) AB674349 AB674386 3 O O 
MEO11 Bulgaria: Central Balkan AB674350  AB674387 1 O 
E. oeme 
OE000 Swiss: Spillmatten, Wallis AB674351 AB674388 1 
OE010 Bulgaria: Berkoviska Planina, West Balkan AB674352 AB674389 2 O 
OE011 Bulgaria: Berkoviska Planina, West Balkan AB674353 AB674390 ] 〇 
OE012 Bulgaria: Pomporovo Resort, Rodopi AB674354 AB674391 3 O O 
OE013 Bulgaria: Pomporovo Resort, Rodopi AB674355 AB674392 1 O 
E. ottomana 
OT000 Italy: Monte Baldo AB674356 AB674393 2 
OT010 Montenegro: Durmitor AB674357 AB674394 4 Q: © 
OTO11 Bulgaria: Orelyak, Pirin AB674358 AB674395 2 O 
OT012 Bulgaria: Granchar-Chalet Pass, Rila AB674359 AB674396 2 〇 
OT013 Bulgaria: Pomporovo Resort, Rodopi AB674360 AB674397 3 O 
OT014 Bulgaria: Pomporovo Resort, Rodopi AB674361 AB674398 1 OQO 
E. pandrose 
PA000 Italy: Piemonte AB674362 AB674399 1 
PA010 Bulgaria: , Borovets, Rila AB674363 AB674400 4 O 
E. rhodopensis 
RH010 Bulgaria: Borovets, Rila AB674364 AB674401 8 C) QU Q 
RHO11 Bulgaria: Granchar-Chalet Pass, Rila AB674365 AB674402 1 © 
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Table 2. Genetic distances as shown by the numbers of base substitutions. (a) Alps and 
Balkan Peninsula (b) Honshu (Japanese islands) and Continent. 


(a) 


Spesies 


Alps-Balkan 





Erebia pamdrose 6 
Erebia oeme 
Erebia alberganus 
Erebia ligea 
Erebia ottomana 
Erebia medusa 
Erebia euryale 


— NNNNA 


Erebia aethiops 


(b) 





Species Continent-Honshu Locality in Continent 
Carterocephalus palaemon 7 Russia: Khabarovsk 
Anthocharis cardamines 6 Russia: Komsomol'sk na Amure 
Colias palaeno 5 Russia: Yakutia, Ust -Nera 


DNA (ca 40 ng) with reaction constitution as recommended 
by Applied Biosystems. Thermal cycling profile was one 
step at 95°C for 10 min; 25 cycles at 95°C for 30 sec, 45°C 
for 30 sec and 72°C for 1 min 30 sec; one extension step 
at 72°C for 8 min 30 sec. PCR products were processed 
after Exo-SAP treatment. The nucleotide sequences of the 
amplified DNA fragments were determined by the direct 
sequencing method with Big Dye Terminator v1.1 Cycle 
Sequence Kit. As primers for PCR and direct sequencing, 
we used V1 (5'-CCT GTT TCT GCT TTA GTT CA-3’), 
C2 (5’-ATC YTT WGA ATA AAA YCC AGC-3^(Yagi et 
al., 1999) for NDS, and mtD6 (C1-J-1718; 5’- GGA GGA 
TTT GGA AAT TGA TTA GTT CC-3^, Nancy (CI-N- 
2191; 5'-CCC GGT AAA ATT AAA ATA TAA ACT 
TC-3^ (Simon et al., 1994), or LCO (5’- GGT CAA CAA 
ATC ATA AAG ATA TTG G -3’), HCO (5’- TAA ACT 
TCA GGG TGA CCA AAA AAT CA -3’) (Folmer et al., 
1994) for CO1. Sequences were edited and aligned using 
SeqScape V. 2.11 (ABI, Weiterstadt, Germany). We 
conducted fragment comparison of 432 bp for ND5, 510 
bp for CO1, and 942 bp in total. A haplotype analysis for 
the gene dataset and construction of a haplotype network 
were conducted using TCS v. 1.21 (Clement et al., 2000). 


Results 


The ND5 and CO/ haplotypes detected are shown in Table 
1, and haplotype network trees obtained from the sequences 
of each species are shown in Fig. 4. E. ottomana was the 
only species having a variable haplotype in each mountain 
range. Other species were roughly divided into two lineages, 
the northern (Balkan Mountains) and the southern (Pirin, 
Rila, Rodopi). Only E. rhodopensis exhibited a very few 


haplotypes despite the considerable sample numbers 
analyzed, with no significant differences in haplotype 
diversity among the other species. 


The following three patterns were noticed by comparing 
the genetic distance among populations within Bulgaria, 
with the genetic distance between the populations in 
Bulgaria and those in the central Europe. 


(1) Genetic distance is generally proportional to 
geographical distance 


(2) Genetic distance is markedly large, without any clear 
correlation with geographical distance 


(3) Genetic distance is very large between the Bulgarian 
and central European population 


Discussion 


1 Three patterns observed in genetic distance between 
Bulgarian and central European populations 


(1) Pattern of genetic distance generally proportional to 
geographic distance 


This is a pattern in which a distance-based isolation 
mechanism (island model) is presumed to be functioning, 
with genetic exchange observed among neighboring 
populations but not among distant populations. In this 
group, Erebia aethiops, E. euryale and E. medusa are 
included. These species are all widely distributed and 
inhabit grasslands in coniferous forest zones. The elevation 
of habitats confirmed by the present authors' field survey 
were 950-1,000 m (West Balkan) and 1,460 m (Rodopi) 
for E. aethiops and 950-1,000 m (West Balkan) for E. 





NII-Electronic Library Service 





The Lepidopterological Society of Japan 


Phylogeography of the genus Erebia in Bulgaria 3] 


'e。 


a) Erebia aethiops 





c) Erebia medusa 
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e) Erebia oeme 


9 


9) Erebia alberganus 





(i) Erebia rhodopensis 


e 





(b) Erebia euryale 





d) Erebia ligea 







(f) Erebia ottomana 


c0 


h) Erebia pandrose 








Fig. 4. Haplotype networks of mitochondrial ND5 and CO/ sequences. 
The size of circles indicates sample size of each haplotype (small circle, n S 2, large circle, n > 2). One 
node indicates a single base substitution and a small filled circle indicates a missing haplotype. 
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euryale. E. medusa was only found in alpine meadows at 
1,680-2,000 m (Central Balkan) and 2,380-2,400 m (Pirin) 
in July, although it is observed in places at around 500 m 
in its early (May) appearance in Bulgaria (Abadjiev, 1993). 


This pattern is easily understood based on the biological 
characteristics of these species. The distribution ranges 
of these species would be expanded together with the 
expansion of the forest zone in the postglacial period, 
while they were fragmented in the ice age. 


(2) Pattern in which genetic distance is markedly large, 
with no correlation with geographical distance 


This pattern would be caused by two populations forming 
a hybrid zone or facing each other in a close distance 
resulting from the subsequent range expansions after the 
reproductive isolation in multiple refugia. Some unknown 
phylogeographical factors were suggested. E. oeme and 
E. ligea are included. 


E. ligea inhabits similar environments to those of E. 
aethiops and E. euryale of the above mentioned group. 
While it often is observed sympatrically with these species, 
there is an extremely large genetic distance between 
populations in northern and southern Bulgaria. In addition 
the Balkan Mountains haplotype is genetically very closely 
related to Switzerland populations, and the south Bulgarian 
populations to geographically distant central Asian 
(Mongolian) populations. E. ligea is widely distributed 
from Europe along the mid-latitude regions of the Eurasian 
continent to the Russian Far East and the Japanese islands. 
The implication is that a ancestral panmictic population 
at an early glacial period divided into two separate 
populations of the eastern lineage and western lineage as 
a result of isolation in more than two refugia during the 
following warm interglacial period. It is presumed that 
these two lineages expanded their distribution areas 
followed by confronting each other in the central Balkan 
Peninsula. As regards populations in the Japanese islands, 
about 15 base substitutions are recognized compared with 
the Mongolian population. More samples from various 
points covering the whole distribution ranges are necessary 
to clarify its phylogeographic history. 


E. oeme, endemic to the European regions, is limited to 
grasslands in coniferous forests in the central European 
Alps. In Bulgaria many populations are also widely 
distributed in the grasslands above the timberline. The 
distribution range is spread from the Pyrenees in the west 
to the Carpathians and Balkan Peninsula in the east, and 
many subspecies and forms are recognized. Although the 
populations in the Balkan Peninsula together with those 
from eastern Austria are classified as the same subspecies 
spodia, having a larger wing span and well developed 


orange-red markings, there is a great genetic distance 
between the populations in northern and southern Bulgaria, 
with no correlation to the geographic distances and 
similarity of morphological characters. Two lineages in 
northern and southern Bulgaria might have been isolated 
in separate refugia prior to the last ice age. 


(3) Pattern with larger genetic distance from central 
European populations 


This group includes E. pandrose and E. alberganus. These 
species have large genetic distances from the central 
European populations, but a comparison was not possible 
between the north and south Bulgarian populations in this 
study. The distribution ranges of each species might have 
been fragmented prior to the last ice age, and apparently 
there might have been no genetic exchanges in or after 
the following cool period up to the present time. Both 
species inhabit grasslands near the timberline, and E. 
alberganus also inhabits grasslands in coniferous forests 
at a lower elevation, but the distribution ranges of these 
species are sporadic. Species having a similar distribution 
pattern have presumably remained genetically isolated 
from prior to the last ice age until today. 


2 Species endemic to the Balkan Peninsula 


Species endemic to the Balkan Peninsula include E. 
rhodopensis, E. orientalis and E. melas, and among them 
only E. rhodopensis and E. orientalis were analyzed in 
this study. The Balkan Peninsula serves as the main 
distribution area for E. ottomana, which was also treated 
as belonging to this group. 


Closely related to E. aetheopella in central Europe, E. 
rhodopensis is mainly distributed in the shrub zone near 
the timberline and also inhabits in a stony or rocky grassland 
above 2,650 m elevation in Bezbog Chalet (Pirin Mts) and 
Grqnchar Chalet (Rila Mts). Genetic variation was not 
detected in this study. 


E. ottomana belongs to the tyndarus group which is well 
diversified in Europe. Its main distribution area is the 
Balkan Peninsula and is also locally distributed in 
southeastern France and northern Italy. The widely 
distributed haplotype was found from the populations of 
Balkan Mountains in the north and Rila Mts in the south, 
and also detected from the population in Montenegro. The 
star-shaped network tree (Fig. 4) contains this widely 
distributed haplotype at the center and a single-base 
substitution haplotype distributed in the Rodopi Mts and 
two-base substitution haplotype distributed in the Pirin 
Mts and central Europe (Monte Baldo, Italy) are derived 
from the main one. Within the tyndarus group, this species 
belongs to a relatively ancient lineage with E. iranica and 
E. graucasica distributed in the Near East, and seems to 
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have expanded its distribution from the Balkan Peninsula 
to central Europe (northern Italy and southeastern France). 


3 History of genetic differentiation 


Here we dare to discuss the relatedness between the 
geographical age when populations of some Erebia species 
in the Balkan Peninsula were reproductively isolated from 
the central European populations and that when the Japanese 
alpine butterflies were genetically isolated from populations 
on the continent. The Balkan Peninsula has been a part of 
the European continent all through the Quaternary period 
(Hewitt, 2000). Honshu, the main island of the Japanese 
islands, however, has been separated from Hokkaido to 
the north by the Tsugaru Strait and from the Korean 
Peninsula to the southwest by the Tsushima Strait. These 
two straits are deep. Honshu separated from Hokkaido 
and the Korean Peninsula and became an island during 
the interglacial age after Riss glaciation about 150,000 
years ago. Itis believed that Honshu has not been connected 
with the continent since then (Oshima, 2000). In other 
words, alpine organisms which now show isolated 
distribution in the alpine belt of Honshu had expanded 
their distribution from the continent prior to about 150,000 
years ago, and they are the species which were affected 
by the climatic changes since then during the interglacial 
period and the last glacial period. At least two periods, 
when glaciers developed, existed in the high altitude areas 
of central Honshu in the Japanese islands, and even now 
there are mountainous permafrosts in some limited areas 
in the high mountainous areas in the Japanese islands 
(Iwata, 2003). 


The same base sequences of ND5 and CO1 regions used 
in this study were also used in the previous studies on the 
genetic differences of the Japanese alpine butterflies from 
the populations on the continent (Nakatani et al., 2005, 
2007a, 2007b), and it allows us to compare the geographical 
age when the Balkan populations were reproductively 
isolated with that when the Japanese alpine butterflies 
were isolated by using the base substitution number. 


The base substitution numbers detected are summarized 
in Table 2. The base substitution numbers between the 
central European and Bulgarian populations are 4-6 for 
the sporadically distributed species: Erebia pandrose and 
alberganus (Table 2a). For the alpine butterflies in Honshu, 
Anthocharis cardamines, Colias palaeno, and Caterocephallus 
palaemon, the base substitution numbers between 
continental and Honshu populations are 5-7 (Table 2b). 
These data suggest that the Bulgarian populations of the 
sporadically distributed Erebia species group were 
reproductively isolated from the central European 
populations during the cooler period prior to the last glacial 
age. This implies that not all northern butterflies were able 


to expand their distribution even in the cooler periods on 
the continent. It also suggests that similar climatic 
environments would have caused such a disjunct distribution 
for northern butterfly populations during almost the same 
period both in Europe and in eastern Asia. 


On the other hand, the base substitution numbers between 
the central European and Bulgarian populations are 1-2 
for the widely distributed species E. aethiops, E. euryale 
and E. medusa. This indicates that genetic exchanges 
between the central European and Bulgarian populations 
of these species had continued up until the last glacial 
period. In addition, genetic exchanges continued until very 
recently between the northern and southern populations 
in Bulgaria for each of these three widely distributed 
species and in addition for E. ottomana originating in the 
Balkan Peninsula and for E. rhodopensis endemic to the 
Balkan Peninsula. 


The vast phylogenetic divergence of E. ligea in the Balkan 
Peninsula is extremely interesting. It is important to clarify 
the genetic structure of this species by analyzing more 
samples from various places covering the entire main 
distributional areas in a future study. 
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ブル ガリ ア に お ける ベニ ヒ カ ゲ 属 蝶 類 の 系 統 地 理 
(中 谷 貴 書 ・ 宇 佐美 真一 ・ 伊 療 建 夫 ) 


バル カン 半島 は , イベ リア 半島 , イタ リア 半島 と 共に 氷河 
時 代 に は , 南方 系 生物 に と っ て レフ ュ ジ ア と し て 機能 し て 
いた た め , ヨー ロッ パ 産 生物 多様 性 の ホッ トス ポッ ト と し 
て 知ら れる . バル カン 半島 と は 一 般 に ドナ ウ 川 より 南 を 指 
す . 西側 に は Dinaric Alps が 北西 か ら 南 東方 向 に 連なり , 最 
高峰 は Prokletije (2,692 m) で ヨー ロッ パ で 五 指 に 入る 大 山 











NII-Electronic Library Service 


The Lepidopterological Society of Japan 





Phylogeography of the genus Erebia in Bulgaria 35 


脈 で ある . その 南端 で は 東西 に 連なる Balkan Mts (最高 峰 
は Betov, 2,376 m) と , 南北 に 連なる Pindos Mts (最高 峰 は 
Katafidi, 2,393 m) に 分 か れる . 蝶 類 は ブル ガリ ア 214 種 , 
ギリ シア 232 種 に 達し , 地域 の 面積 か ら す る と 上 虹 類 の 種 数 
で る も 多様 性 が 高い こと が 伺わ れる . 


サン プル と し て , ブル ガリ ア の 北部 (West and Central Balkan 
Mts) お よび 南部 (Rodopi, Pirin, Rila Mts) で 採集 し た Erebia 
属 10 種 65 個体 , お よび スイ ス , イタ リア な ど 中 欧 . お よび 
中 央 アジ ア 産 の 10 種 11 個体 を 用 いた . ミト コン ドリ ア の 
遺伝 子 領域 と し て は , 種 内 変異 の 解析 に よく 使わ れる NDS 
の 一 部 (432 bp) と CO7 領域 の 一 部 (510 bp), 合計 942 bp 
を 使用 し , 得 ら れ た 塩基 配列 を 比較 し て ハプロタイプ を 決 
定 し た . 


ブル ガリ ア 国 内 産 と 中 欧 産 個体 群 間 の 遺伝 距離 に み ら れ る 
パタ ー ン 
ブル ガリ ア 国 内 に お ける ハプロタイプ の 分 化 と 中 欧 産 と の 
遺伝 距離 を 比較 する と , 3 つの パタ ー ン に 分 類する こと が 
で きる . 


(1) 遺伝 距離 が ほ ば 地理 的 距離 に 比例 し て いる パタ ー ン 


近隣 個体 群 間 で は 遺伝 的 交流 が 認め られ る が , 遠 距離 の 個 
体 群 間 で は 認め られ な い , いわ ゆる 距離 に よる 隔離 機構 
( 島 モ デル ) が 働い て いる と 想定 され る パタ ー ン で ある . 
この グル ー プ に 属す る 種 と し て は , E. aethiops, E. euryale, 
E. medusa EA. いずれ も 針葉樹 林 淀 の 草原 に 生息 する 広 
域 分 布 種 で . その 生態 的 特性 か ら も 推定 され る パタ ー ン で 
あっ て , 氷河 期 に 分 布 が 分 断 さ れ た も の が , 温暖 期 に お け 
る 森林 淀 の 拡大 に 伴っ て 分 布 を 広げ た と 推定 され る . 


(2) 個体 群 間 の 地理 的 距離 と は 相関 の な い 遺 伝 距離 が 認め 
られ る パタ ー ン 


過去 に 複数 の レフ ュ ジ ア に 生殖 隔離 され る こと で 固定 され 
た 複数 の 系 統 が , その 後に 分 布 を 拡大 し て 混 生 地 を 生じ た 
り , ある い は 地理 的 に 至近 距離 で 対 財 し て いる 分 布 型 で あ 
り , 何ら か の 系 統 地 理 的 要因 の ある こと が 示唆 され る . E. 
oeme, E. ligea が 属す る . 


E. jeeg は 前 の グル ー プ に 属す る aethiops, E. euryale と 類 
似 の 環境 に 生息 し , し ば し ば これ ら の 種 と 同 所 的 に 見 られ 
る 場合 が 多い に も か か わら ず , ブル ガリ ア 北 部 の バル カン 
山脈 産 と 南部 の 山系 産 の 個体 群 間 の 遺伝 距離 は きわ め て 
大 きい . さら に バル カン 山脈 産 の ハプロタイプ は スイ ス 産 
個体 群 と ,. また ブル ガリ ア 南 部 山系 産 個 体 群 は 地理 的 に 遠 
< く 離れ た 中 央 アジ ア (モン ゴル ) 産 の 個体 群 と 遺伝 有 的 に 極 
め て 近 縁 で ある こと が 判明 し た . E. ligealk 8 7 8 v» 2S5 
ュ ユーラシア 大 陸 の 中 緯度 地域 に 沿っ て ロシア 極東 , 日 本 列 
島 に まで 広域 分 布 す る 種 で ある が , 東西 2 系 統 に 分 断 さ れ 
た 個体 群 が . その 後に 分 布 を 拡大 し バルカン 半島 の 中 部 で 
両 系 統 が 対 財 し て いる と 想定 され る . 今後 大 陸 各 地 に お い 
て サン プル 地点 を 増やす こと で 系 統 地理 的 な 全貌 を 解明 す 
る 必要 が ある . 


E. oeme は ヨー ロッ パ 地 域 特産 で , 中 欧 の アル プス で は 針 
葉 樹林 内 の 草原 に 生息 し 分 布 は 局 限 さ れる の に 対し て , T 
ル ガ リ ア で は 森林 限界 より 上 の 草原 に も 広く 分 布 し 個体 数 
も 多い . 西 は ピレネー 山脈 か ら 東 は カル パチ ア 山 脈 ・ バ ル 
カン 半島 に 分 布 し . いく つか の 亜 種 また は フォ ー ム に 区 分 
され る . バル カン 半島 産 は 東 オ ー ス トリ ア 産 と 共に 大 型 で , 
斑紋 が よく 発達 する 特徴 を 持ち 亜 種 ゅ oZjg と され る が , At 
部 ブル ガリ ア 産 と 南部 ブル ガリ ア 産 個体 群 の 間 の 遺伝 距離 
は 大 きく , 遺伝 距離 と 形態 的 特徴 の 類似 性 に は 相関 性 が な 
い . 両者 の 個体 群 は . その 遺伝 距離 の 大 きる さか ら み る と 最 
終 氷河 期 以前 か ら 別 々 の レ フュ ジア に 隔離 分 布 し て いた こ 
と が 示唆 され る . 


(3) 中 欧 産 個体 閣 と の 遺伝 的 距離 が 大 きい パタ ー ン 


中 欧 産 個体 郡 と の 遺伝 距離 が 大 きい グル ー プ で ある . 最終 
水 期 よ り も 前 に 分 布 が 分 断 さ れ , その 後 の 寒冷 期 に も 遺伝 
的 交流 が な いま ま 現 在 に 至っ て いる と 想定 され る . ち 
pandrose, E. alberganus が 属す る . E. pandrose は 森林 限界 付 
近 の 草原 に 生息 する が , E. alberganus は さら に 低 標 高 の 針 
葉 樹 林内 の 草原 に も 生息 する . スイ ス で は 高山 草原 で も 草 
丈 の 高い 草原 に の み 見 られ , 食 草 が 限定 され て いる の か も 
し れ な い . この よう な 離散 的 分 布 域 を も つ 北 方 系 蝶 類 の 一 
部 の 種 は , 最終 氷 期 より も 前 に 分 布 が 分 断 さ れ . TORO 
寒冷 期 に も 分 布 を 拡大 する こと が で き な か っ た と 推定 され 
る . 


バル カン 半島 固有 種 


バル カン 半島 に は 万 rhodopensis, E. orientalis, E. melas の 
固有 種 が 知ら れ て いる が , これ ら の 内 E. rhodopensis, E. 
orientalis に つい て 解析 で きた . バル カン 半島 が 主たる 分 布 
地 の otomana も この グル ー プ に 属す る も の と し て 扱う . 


E. rhodopensis は 中 欧 に 分 布 する 万 aetheopella E Iii% T H 
り , GR UR FL PS IED HEA GH IC PATO HUD 2345 V, Bezbog 
Chalet (Pirin Mts) や Grqnchar Chalet (Rila Mts) で は 標高 
2650 m EA LO BMT CSAS. 調査 し た 範囲 で は 遺伝 
的 変異 は 認め られ な か っ た . 


E. ottomana は 中 欧 で 多様 な 種 分 化 を 遂げ た py グルー 
プ に 属し , フラ ンス 南東 部 と イタ リア 北部 に も 分 布 す る が , 
バル カン 半島 が 主たる 分 布 城 で ある . 広域 分 布 型 の ハブ プロ 
タイ プ が 北部 バル カン 山脈 と 南部 Rila 山系 か ら 見 出さ れ , 
同じ 型 は モン テ ネ グ 口 産 個体 群 か ら も 検出 され た . ネッ ト 
ワー ク 樹 に 示す よう に , この 広域 分 布 弄 ハプロタイプ を 中 
心 に し て , 1 塩基 置換 の ハプロタイプ が Rodopi Mts に , 2 Hit 
基 置 換 の タイ プ が Pirin Mts と 中 欧 (イタ リア , Monte Baldo) 
に 分 布 す る スタ ー 型 の 構造 を 示す . 本 種 は rndarus グル ー 
プ の 中 で は 近東 に 分 布 す る iranica, E. graucasica と 共に 
古い 系 統 に 属す る 種 で あっ て , バル カン 半島 か ら 中 欧 ( イ 
タリ ア 北 部 , フラ ンス 南東 部 ) へ 進出 し た も の で あろ う . 


E. orientalis は Rila Mts の みか ら 得 られ た た め . ブル ガリ ア 
南北 の 個体 群 間 の 比較 は で き な か っ た . 
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遺伝 的 分 化 の 歴史 これ ら の 種 の 中 欧 と バル カン 半島 の 個体 群 が , 最終 氷 期 を 


日 本 産 高山 農 の 各種 に つて は , 筆者 ら に よる 大 陸 産 の 同種 
と 遺伝 的 差異 を 比較 検討 し た 研究 ., お よび 日 本 列島 内 で の 
ハプロタイプ を 解析 し た 研究 が ある . 本 州 は 最終 氷 期 に は 
大 陸 や 北海 道 と 陸続き に な ら な か っ た た め , 本 州 に 遺存 分 
布 する 高山 蝶 は 最終 氷 期 よ り 前 の 寒冷 期 に 大 陸 か ら 渡 来 
し た も の と 考え られ る . し た が っ て 本 州 の 例 と 比較 する こ 
EC, バル カン 半島 で の 生 猪 隔離 が 始ま っ た 年 代 を 推定 す 
る こと が で きる . 


分 析 し た 遺伝 領域 MDs5 と CO7 の 塩基 置換 数 を まとめ る と 
Table 2 の よう に な る . 中 欧 と ブル ガリ ア 産 個体 群 間 の 塩基 
置換 数 は , 分 布 が 離散 的 な pandrose と E. alberganus は 4 
~6 で ある (Table 2a). ~ER EBED EYY VS 
チョ ウ , ミヤ マモ ン キ チョ ウ , タカ ネ キ マダ ラ セ セ リ は , 
大 陸 産 と 本 州 産 個 体 群 の 塩基 置換 数 が $ 一 7 で ある (Table 
2b). これら の デー タ か ら , HEB AO Erebia RESH 
欧 と プル ガリ ア 間 で 生殖 隔離 し た 時 期 は , 最終 水 期 よ り 前 
の 寒 傘 期 で ある と 示唆 きれ る . すなわち ヨー ロッ パ と アジ 
ア 東 部 に お いて , ほぼ 同時 期 に 北方 系 蝶 類 の 個体 群 が 隔 離 
分 布 す る よう な 気候 環境 で あっ た こと が 示唆 され る . また 


通じ て 遺伝 的 交流 が な か っ た こと は , 寒冷 期 と は いえ ぇ 北方 
系 蝶 上 類 の すべ て が 分 布 域 を 拡大 で きた わけ で は な いこ と を 
示唆 し て いる . 


同じ く Table 2a に よる と , 中 欧 と ブル ガリ ア 産 個体 群 間 の 
塩基 置換 数 は , 広域 分 布 型 の Erebia aethiops, E. euryale, E. 
medusa で は 1 一 2 で あり , 最終 水 期 以降 も 遺伝 的 交流 が あ 
っ た こと が 示唆 され る . また ブル ガリ ア の 南北 個体 群 間 で 
は , 広域 分 布 型 の Erebia aethiops, E. euryale, E. medusa, お 
よび バル カン 半島 が 発祥 の 地 と 推定 され る 万 ottomana, バ 
ルカ ン 半 島 固 有 種 の ち rhodopensis V2 に お いて は , E. aethiops 
以外 は 広域 分 布 型 の ハプロタイプ が 南北 の 山 域 に 共通 し て 
生息 し て お り , ご く 最 近 ま で 遺伝 的 交流 の あっ た こと が 示 
唆 され る . 

バル カン 半島 で 系 統 分 岐 の 推定 され た 互 ligea は 系 統 地理 
的 に 極め て 興味 深い 事象 で あっ て , 今後 さら に サン プリ ン 
グ 地 点 を 増やし て , 主要 な 分 布 域 を カバ ー す る 集団 の 遺伝 
的 構造 を 明らか に する 必要 が ある . 
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